Momentum-resolved study of an array of ID strongly phase-fluctuating Bose gases 
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We investigate the coherence properties of an array of one-dimensional Bose gases with short- 
scale phase fluctuations. The momentum distribution is measured using Bragg spectroscopy and 
an effective coherence length of the whole ensemble is defined. In addition, we propose and demon- 
strate that time-of-flight absorption imaging can be used as a simple probe to directly measure the 
coherence-length of ID gases in the regime where phase-fluctuations are strong. This method is 
suitable for future studies such as investigating the effect of disorder on the phase coherence. 
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Physics of one-dimensional (ID) systems attracts a 
great interest both on the theoretical and the experimen- 
tal side. Recent progress in nanotechnology allowed to 
implement ID systems in a variety of fields, from inor- 
ganic and organic superconductors Q , carbon nanotubes 
and nanowires 1 to spin chains and ladders |, as well 
as cold atomic systems [fl 0] • All these systems belong to 
the universality class of interacting quantum fluids known 
as Luttinger liquids [§], whose properties strongly dif- 
fer from their 3D counterparts. For instance, quantum 
and thermal fluctuations are strongly enhanced by the 
reduced dimensionality, their knowledge giving access to 
key quantities characterizing the system [9] . Their pres- 
ence can drastically alter the properties of the systems, 
such as in superconductive disordered nanowires where 
they can lead to the formation of phase-slip centers [l(| ■ 

In the context of cold atoms, both phase and density 
fluctuations of ID systems have been studied in the last 



years [llHl4j. In particular, phase coherence has been 
investigated by monitoring interference between two dif- 
ferent ID gases [ijj and by observing density modula- 
tions [3] or the response to light scattering [3] in elon- 
gated 3D quasi-condensates. Nevertheless, in all these re- 
alizations transverse trapping frequencies hardly exceed 
a few kHz and are typically of the order of the chemical 
potential and the temperature. Reaching the regime of 
strongly interacting ID systems would further enhance 
the presence of quantum and thermal fluctuations. This 
can be achieved when atoms are loaded in 2D optical lat- 
tices allowing for much stronger transverse confinements 
d, 0] • Yet in the latter case one obtains a larg e num- 
ber of ID tubes for which techniques like 12], [l3| cannot 
be implemented to study the phase coherence properties 
since tube-averaging washes out the response signal. 

In this work we investigate the axial coherence prop- 
erties of an array of strongly phase-fluctuating ID Bose 
gases and we suggest time-of-flight (TOF) imaging as a 
probe of the coherence length. In our case, thermally- 
induced phase-fluctuations dominate and drastically re- 
duce the coherence length (2L^) of the system compared 
to a 3D Bose- Einstein condensate (BEC) [ljj . We use 



Bragg spectroscopy with large momentum-transfer [161 ] 
to measure the momentum distribution and directly eval- 
uate Lff, [3]. In addition, we verify that direct mapping 
of the momentum distribution into coordinates-space via 
absorbtion imaging after TOF is an effective probe of the 
phase-fluctuations. We demonstrate that in our range of 
parameters these two techniques yield the same results. 
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FIG. 1: (a) Schematic view of the two-dimensional array of 
one-dimensional gases, (b) Bragg spectrum of an array of 
strongly correlated ID gases produced by lattices with am- 
plitude s — 50. Red curve is a fitting of the function vf(y), 
where f(y) is a Lorentzian. (c) Half width at half maximum 
of the Bragg resonances as a function of the amplitude s of 
the lattices, for two different value of momentum [blue squares 
(qi = 16.0(2)/im _1 ) and black dots (q 2 = 7.3(2)/zm -1 )]. 

Our system, sketched in Fig. [Ha), consists of about 
2x 10 3 ID atomic micro-tubes. Each gas has typical total 
size ^30 /im x 0.05 /im and linear density n\u ~5/im -1 . 
To arrange atoms in this configuration, we confine a 3D 
BEC of 87 Rb in a pair of orthogonal red-detuned optical 
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lattices. We study different configurations by tuning the 
amplitude V of the 2D optical lattice (s — V/Er ranging 
from to 56, where Er = h? /2m\ 2 L , h being the Planck's 
constant, m the atomic mass and A^ = 830nm the lattice 
wave-length). The stronger the optical confinement, the 
more anisotropic is the trap experienced by each ID gas 
(the aspect ratio A = cJx/^ih namely the ratio of the ra- 
dial harmonic trapping frequency to the axial one, ranges 
from 787 to 880 for 5 < s < 56). For all the amplitudes 
of the transverse lattice we explore, each gas has a fully 
ID character, i.e. both chemical potential and tempera- 
ture are about one order of magnitude smaller than the 
frequency of the transverse harmonic oscillator. The cru- 
cial quantity to describe the regime of the ID gas is the 
parameter 7 = rngm /^ 2 ^id, that is the ratio of interac- 
tion energy to the kinetic energy necessary to correlate 
particles at distance 1/nuj, gio being the interatomic 
coupling in the ID gas. In our case 7 ~ 0. 3-1-0. 6 so 
that inter-particle correlations are stronger than in the 
mean- field regime (see for example Q). 

We first investigate the effect of the phase-fluctuations 
via Bragg spectroscopy. In brief, the lattice gas is 
diffracted from a moving lattice created by two simul- 
taneous off-resonant light pulses (Bragg beams) with a 
relative angle 9, detuned from each other by a tunable 
frequency-difference w/(2ir). This perturbs the system 
activating excitations with energy fnu and momentum 
fki B , the modulus of which depends on 9 [l6j]. The ge- 
ometry of the Bragg beams is chosen to align hq B to the 
axis of the ID tubes. After the excitation, the lattice 
amplitude is turned down in 15ms to a low value (s = 5) 
where the different tubes are no more independent allow- 
ing the system to re-thermalize via atom-atom collisions. 
After 5ms both optical and magnetic trap are simulta- 
neously switched off and the system is observed after a 
time of flight (TOF) t T oF = 21ms. The physical ob- 
servable is the increase of the size of the central peak of 
the atomic cloud. More detailed comments on the exper- 
imental procedure have been reported in Ref.[l7|]. Our 
technique allows to measure the energy AE transferred 
to the system [18] , which depends on the imaginary part 
of its polarizability x"( w )j apart from the characteris- 
tics of the perturbing potential (amplitude Vb and time 
duration At) [l9j]: AE cx V|wx"(w)At. The polarizabil- 
ity can be expressed in terms of the dynamical structure 
factor of the system = xS(u},q)(l - e -^/( fc « T )). 

Two different geometrical configurations of the Bragg 
beams have been used to vary the transferred momen- 
tum along the axis of the gases. Counter-propagating 
beams along the axis of the atomic tubes yield q% = 
16.0(2)^m , a small-angle configuration gives q 2 = 
7.3(2)^m -1 . In both the cases we assume the excita- 
tion to be in the Doppler regime [20| . where S(a, co) 
is reduced to the momentum distribution n(q) flit Il6|. 
In this regime the spectral half-width at half-maximum 
(HWHM) can be related to the momentum width HAq 



through the relation Av — i^—hAq, which is linear in the 
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wave- vector qs of the excitation [21] . In the experiment, 
the ratio of the HWHMs of the response of identical ar- 
rays of ID gases to the two different excitations qi and 
q 2 is consistent with q\/q% = (2.16 ± 0.06), as expected 
(linear fitting of the experimental data on Figfljc) allows 
for defining a mean ratio Av\/Av^ = (2.7 ± 0.8)). 

In our range of 7 w 0.3 4- 0.6, interactions are beyond 
the mean-field description but still far from the Tonks- 
Girardeau regime. Thus, we expect the interaction- 
induced spatial decay of one-particle correlation function 
to happen on a larger scale than that leaded by phase- 
fluctuations, for typical temperatures in the experiment 
(^100nK) [22j. The one-particle correlation function be- 
ing dominated by the exponential decay due to phase- 
fluctuations, the momentum distribution exhibits a pro- 
file well described by a Lorentzian shape [l5j . Its HWHM 
HAq relies only on the coherence length of the gas, being 



Aq = 



0.635 



(1) 



where = K 2 niD/ '(mfc^T) is the half coherence-length 
(T being the temperature) ; the factor 0.635 take into 
account the ID geometry [231 ] . 
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FIG. 2: L$ of the array of ID gases is reported as a function 
of the amplitude s of the lattices which squeezes the gas in 
ID micro-tubes and for two different values of the momentum 
of the excitation [blue squared (qi) and black dots ((fc)]- 

The description of the problem is complicated by 
the presence of an array of gases with different densi- 
ties (thus, different characteristic L^). In principle one 
should consider the sum of the response of each tube. 
Supposing a mean-field picture [2^ . interactions would 
give a broadening of the width and a shift of the center 
of the Bragg resonance compared to the single-particle 
response hv s , v , — fi 2 q 2 /(2m), both depending on the 
density of each tube [l6|. Yet the global response of 
the system to the Bragg excitation consists in a single 
broad resonance, as depicted in Fig. QJb), the center of 
which is shifted compared to v s . p . (isq > v s .p.)- Its shape 
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is well described by vfiy) — LuS(q,ui) where f(v) is a 
Lorentzian function (see Fig. [TJb)). This suggests that 
thermal broadening of the response of each gas exceeds 
the interaction-induced broadening and masks the rela- 
tive shifts of the resonant frequencies of the tubes. We 
verified numerically that it works for our experimental 
parameters 25| . Therefore we analyze the Bragg spectra 
as being the response of a single ID gas and we define 
accordingly an effective coherence-length of the whole 
system, using the relation in Eq. [T] 

From the fittings of the Bragg spectra, we extract the 
HWHM Aza This quantity is reported in FigQJc) as a 
function of the amplitude s of the optical lattices and 
for the two different wave- vectors of the excitation q\,q2- 
The total number of atoms is kept almost constant in 
both the series of data. From the spectral half-width, we 
extract the half-length L^. As shown in Fig 121 rapidly 
drops down as s increases. In addition, the analysis of 
the spectra for the two momenta q\ and q± reveals consis- 
tent as expected in the Doppler regime. The optical 
confinement makes the aspect ratio of the ID gases grow 
and their density decrease. However, we estimate numer- 
ically |25| that the relative variation of the ID-density in 
the whole range of s is about 10% and does not justify 
the rapid downfall of the coherence length. This suggests 
that major role in determining is played by the finite 
temperature. Note that extracting temperature from the 
measurement of is not straightforward as it requires 
to take into account the inhomogeneity of nm over the 
array; this will be the subject of a future work [25j . 
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FIG. 3: Mean-square value of residuals of a Lorentzian (hallow 
dots) and Gaussian function (solid dots) fitted to the momen- 
tum distribution mapped through TOF. 

To keep the insight up on the system, informa- 
tion on the coherence-length induced by thermal phase- 
fluctuations has been also extracted by direct map- 
ping the momentum-distribution into space-distribution, 
which is measured via absorption imaging of the gas after 
switching-off the trap (below referred to as TOF measure- 
ments). The expansion of the atomic gas from the trap is 
governed by two kinds of kinetic-energy: the one which 



interactions convert into, and the one produced by in- 
trap phase-fluctuations. Due to the strong anisotropy 
of the trap, the interaction-induced expansion mainly 
affects the radial direction [26], whereas the longitudi- 
nal size of the cloud R™q F ^ s n0 * significantly altered 
compared to its in-trap value. At non-zero temperature, 
thermally-induced local phase-gradients produce a veloc- 
ity field given by v,£ = (h/m)W4> [13] j 4> varying signif- 
icantly on a length-scale L<f,. It determines an increase 
R TOF °f the longitudinal size during TOF, which con- 



tributes relevantly if R T of/^TOF > 1> where 



R 



TOF 



TDint 

n TOF 



(2) 



In previous experiments, the product L^Ri^Qp 
amounts typically to 10^mx260^m for elongated 3D 
quasi-condensates 14[ and ~l/xmxl70^m in the case of 
atom-chip experiments 12 1. In both cases, R T q F is neg- 



ligible and the longitudinal length after TOF cannot be 
related to in-trap phase-fluctuations. For our ID lattice 
gases, this quantity is reduced to ~ l/mix27/^m for s = 5 
and it falls even one order of magnitude as the amplitude 
of the optical confinement increases (~ 0.2^mx22/xm for 
s = 56). Thus one expects that the in-trap size still 
dominates for low values of s and the density profile has 
a parabolic shape (possibly smoothed to a Gaussian by 
the finite resolution of the imaging system) , as expected 
in the Thomas-Fermi regime. For high values of s phase- 
fluctuations should enlarge the distribution and the pro- 
file should assume a Lorentzian shape. To confirm this 
behaviour we have analyzed the TOF profiles with both 
Gaussian and Lorentzian functions. From these fittings 
we plot in Fig. |3] the mean squared value of the residuals 
of both fitting functions. As anticipated from the simple 
formula in Eq.[2l Fig. [3] points out the cross-over through 
the two regimes at s ~ 10 4- 15. In particular, for s > 20 
initial trapped phase fluctuations are responsible for the 
Lorentzian shape of the TOF momentum distribution. 

To quantitatively compare the results of Bragg spec- 
troscopy and TOF measurements, we map both energy 
spectra and density profile after TOF into momentum- 
space of the in-trap gas. In the first case, we use the free- 
particle dispersion relation obtaining q = ^n 2 mv / (hqs) — 
qB/%- In the latter case, the calibration of the pixel size in 
momentum-space is obtained by measuring the distance 
between two interference peaks released from the lattices 
at weak amplitude. For an array of strongly-correlated 
ID gases (s = 50 in FigHJa)) the momentum distribu- 
tions measured via Bragg spectroscopy (data points) and 
TOF measurements (continuous line) show an excellent 
agreement. From TOF measurements we extract the co- 
herence length as well. Accordingly to our results from 
Bragg spectroscopy, is observed to reduce as s in- 
creases (Fig. S(b)). In the inset, we compare mea- 
sured in the two ways. The accordance fails in the region 
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of parameters (gray area) where TOF profiles are not 
dominated by phase- fluctuations (see Eq. [2]). 
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FIG. 4: (a) Momentum distribution of an array of ID gases 
in an strongly confining optical lattice (s = 50) measured 
through Bragg spectroscopy (black dots) and direct mapping 
in TOF density profile (red curve), (b) Half coherence length 
Ltj, extracted from TOF measurements is shown as a func- 
tion of s in linear scale. Gray area points out the region 
of parameters where TOF measurements can not be used to 
extract Lj,. Inset: Comparison between L<j, from Bragg mea- 
surements (black dots) and direct mapping (red dots). 

In conclusion, we investigate the coherence properties 
of an array of ID Bose gases by measuring their mo- 
mentum distribution. We observe the latter to have a 
Lorentzian shape as predicted for a single uniform ID 
gas. We define an effective coherence length of the whole 
ensemble and we show its evident reduction as the op- 
tical confinement is increased. Comparing Bragg spec- 
troscopy and direct mapping of momentum into density 
distribution after TOF demonstrates that TOF images 
give access to coherence properties in the presence of 
strong phase-fluctuations. Our work paves the way for 
future studies of the coherence properties in ID geome- 
tries with short coherence lengths, and especially strongly 
interacting disordered systems. Of particular interest 
are strongly interacting disordered systems where the 
role of thermal phase fluctuations in the nature of the 
superconductor-insulator transition is debated (28|. So 
far, only disordered quasi-condensates have been investi- 
gated where it was shown that the contribution of phase 
fluctuations is small (29l |. 
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